Skyrmions, Hadrons and isospin chemical potential 

M. Loewf! and S. Mendizaba|l| 
Facultad de Fisica, 
Pontificia Universidad Catolica de Chile, 
Casilla 306, Santiago 22, Chile. 

J.C. RojaJ] 

Departamento de Fisica, Universidad Catolica del Norte, Casilla 1280, Antofagasta, Chile 

Abstract 

Using the Hamiltonian formulation, in terms of collective variables, we explore the evolution of 
different skyrmionic parameters as function of the isospin chemical potential (/x), such as the energy 
density, the charge density, the isoscalar radius and the isoscalar magnetic radius. We found that 
the radii start to grow very fast for fj, > 140 MeV, suggesting the occurrence of a phase transition. 
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The skyrmionic approach to baryon dynamics is an interesting attempt to discuss the 
occurrence of phase transitions induced by temperature and/or density effects in the hadronic 
sector [Ij. The idea is to study the stabihty properties of such objects by an analysis of their 
mass behavior and the spatial extension, i.e. mean square radius, associated to different 
currents. Recently this idea has been extrapolated also to the analysis of Skyrmions in 
curved spaces 

In a previous article, we discussed the static properties of the Skyrmion solutions in the 
presence of a finite isospin chemical potential /i We were able, among other results, to 
find a critical chemical potential /ic = 222.8 MeV, where the Skyrmion mass vanishes. 

In this article we will extend our analysis by using the Hamiltonian formulation introduced 
by Adkins et al. . Note that the identification of baryons as Skyrmion states emerges only 
in the Hamiltonian formalism. The difference between the static analysis and the present 
one, is that now we are able to look into the energy spectrum of non strange nucleons as 
function of u. In the static approach we cannot distinguish between different nucleon states. 

n 

Using collective coordinates (z = 1..4), following p|, we establish the Lagrangian in 
terms of and dj. Then we infer the Hamiltonian operator which gives us the energy levels 
of the nucleon, the charge densities, the isoscalar mean square radius and the magnetic mean 
square radius. In particular, we emphasize the growing behavior of the mean squared radius, 
associated to different conserved currents, as function of /i, suggesting the occurrence of a 
phase transition. 

The Skyrme lagrangian with isospin chemical potential (/x) is given by 



p2 

+ ^,Tr[{D,U)U^,{DM)U^]\ (1) 
As usual, the isospin chemical potential is introduced through the covariant derivative 

m 

D, = d,-t^[as,U]5,,o. (2) 
In ref yL|], the static solution f/o(r) of (Q) was found using a radial symmetric "Hedgehog" 
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ansatz [3|, finding the minimum of tlie Skyrmion mass witli the parametrization 



Uo = exp [iF{r)aini\ , (3) 

where the cjj's are the Pauh matrices, denotes the spatial unitary vector and F(r) is a 
numerically determined function, that satisfies F{0) = vr and -F(oo) = 0. When the isospin 
chemical potential is turned on, the profile acquires also a dependence on fi. 



F{r) F{r,fi). 

With this profile, the energy density acquires the following form Q 



(4) 



with 



= M^=o 
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where we introduced the dimensionless parameter 
r = eF-j^r and the integrals /2, I a 



(6) 



h = j rf=^fTr[f?,L,]% (7) 

with ^7 = (J3 — Ua^U\ (To is the 2x2 identity matrix and Lj^ = {di^U)W . The variational 
equation for the static Skyrme Lagrangian given in (^Q), allows us to find a numerical solution 
for the profile F{r) ^j. As it was already mentioned, the mass of the Skyrmion vanishes for 
a critical chemical potential. 
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As usual, in order to obtain the hadronic spectra, it is convenient to introduce SU{2) 
collective coordinates A{t) ^, such that 



U = A{t)UoA^it). (8) 

The SU{2) matrix A{t) is parameterized by the Pauli a matrices and the identity o"o, 
according to 

A{t) = ao(t)ao + ia{t) ■ a, (9) 
where the a's obey the constraint 

al{t) + a\t) = l. (10) 

Introducing (jSI) and Q into (H)), a direct (but rather involved) computation leads us to 
the Lagrangian 



L = -M„ + 2A 



do + 



/iOi 



+ 02 + — — + 03 



2 ) 



-M^ + 2A d, + /i 



A,, 



(11) 



where A = {2ti /?,e^F^) A, with 



A = J r"^ sin^ F 
In equation (jllj) we have defined 



1 + 4 ( + 



sin^F 



(12) 



Ai QijClj, 



(13) 



where 
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(14) 
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Explicitly, 



Aq = as, Ai = -02, 

A2 = ai, A3 = -ao- (15) 

Notice that the A has the same functional dependance on F as the equivalent parameter 
defined in ^. In our case F depends also on /z. is the chemical potential dependent 
mass given in (0). 

From the lagrangian (jllll . we get the Hamiltonian 

H = M,- 2A/i2 + ^, (16) 
where the canonical momentum is given through a minimal coupling 

7ri=pi- AXfiAi. (17) 
The Hamiltonian can be expressed in the following way 

2 

H = M^ + ^-fiA,pi, (18) 

OA 

and considering the canonical quantization procedure Pi ^ pi = —iS/Sai, we get 



where I3 is the third component of the isospin operator 



f i f S 5 5 \ 

h = - [ '^OT CtfcT eklmO-l- . (20) 

2 V oak oaQ ba^ j 

Following the usual procedure, we may associate a wave function to the Skyrme Hamil- 
tonian. In order to identify baryons in this model, these wave functions have to be odd, i.e. 
■?/'(A) = —■?/;(— A). In particular, nucleons correspond to linear terms in the a's, whereas the 
quartet of A's are given by cubic terms. 



5 



T 1 1 r 




|x(MeV) 



FIG. 1: The Nucleoli static energy as function of the isospin chemical potential fx. The splitting 
between neutron(solid line) and proton(dotted line) is due to the last term in 

The energy spectra of nucleons as function of /i is shown in figure d We can see that an 
energy sphtting between neutrons and protons is induced. This Hamiltonian remind us the 
Zeeman effect, where the generation of the energy spectra is broken, by an external magnetic 
field. In our case the isospin chemical potential plays the same role. 

Here we are interested to establish the relation between the relevant physical parameters 
and the isospin chemical potential, for that purpose, we will consider the Baryonic, vector 
and axial currents, exploring the behavior of the effective radii associated to them. 

The Skyrmion model allows the existence of different conserved currents and their re- 

n 

spective charges [8||. Using different charge densities we may define several effective radii. 
The evolution of those radii as function of chemical potential provides information about 
the critical behavior close to the phase transition, where the Skyrmion is no longer stable. 
Let us first start with the topological baryonic current 

= [iU^dM){U^d^U){U^dpU)] . (21) 

The baryonic charge density for the Skyrmion is given by 

PS = 47ir^B°{r) = -- sin^ F{r)F'{r). (22) 

TT 

Obviously, drpB = 1, independently of the shape of the skyrmionic profile. 
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FIG. 2: The isoscalar mean square radius as function of ^. 
The isoscalar mean square radius is defined by 



oo 



{r')i=o= / drr^pB. (23) 





This radius seems to be quite stable up to the value of ~ 120 MeV, starting then to 
grow dramatically. Although we do not have a formal proof that this radius diverges at a 
certain critical fi = fic, the numerical evidence supports such claim, as it is shown in figure El 
Divergent behavior for several radii, associated to different currents, has also been observed 
in different hadronic effective couplings as function of temperature in the frame of thermal 
QCD sum rules 0. 

The same kind of behavior is found for the mean square radius that emerges from the 
isoscalar magnetic density 

, r^F'sin^F 

which is plotted in figure El 

This growing behavior of the mean square radius has its counterpart in the electric charge 
distribution. It turns out that the proton charge distribution becomes broader for higher 
chemical potentials. 

Following the usual convention, we introduce the isoscalar and isovector magnetic mo- 
ments 
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FIG. 3: The isoscalar magnetic mean square radius as function of /j,. This figure suggests the same 
divergent behavior as the isoscalar mean square radius with the same critical chemical potential 
(figure HI). 



fIi=o = ^ y ^ X B£x, (25) 
fli=i = IJrx V^d^x, (26) 
where B is the vector part of (1211) and V is the Noether current associated to the vector 

n 

charge. Following |3|], we consider 



2 A 

(/^/=i)3 = ^vr-^. (28) 
From the definition of the Bohr magneton for nucleons 



l^={lL)^^ (29) 



AM. 

and making the identification gi=Q = gp + g-n and gi=i = gp — gni it is possible to obtain the 
magnetic moments for the proton and neutron; /ip = gp/2 and /i„ = 5'n/2. The behavior of 
both magnetic moments is presented in figure El Besides, figure El shows the ratio | /ip/ /i„ | , 
the figure suggests that such ratio goes to one for ii /Xc- 
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FIG. 4: The magnetic moments fip and ^in for the proton, (sohd hne) and the neutron (dotted 
hne) respectively, as function of /i. 
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FIG. 5: The ratio j | as function of /i. 

As a conclusion, we would like to remark the divergent behavior for the isoscalar and mean 
square magnetic radii as function of the chemical potential. This suggests the occurrence 
of a phase transition. In fact in the QCD sum rules approach, the radii are phenomeno- 
logical order parameters for thermal deconfinement Ql . In this case, however, this behavior 
is induced by the isospin chemical potential. Finally, it is interesting to notice that the 
Hamiltonian approach allows us to distinguish between the mass evolution of neutrons and 
protons. 
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